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Abstract 

In the Minimal Super symmetric Standard Model with complex parameters 
(cMSSM) we perform a complete one-loop calculation for the Higgs boson masses 
(including the momentum dependence) and the mixing angles. These corrections are 
obtained in the Feynman-diagrammatic approach using the on-shell renormalization 
scheme. The impact of the newly evaluated corrections is analyzed numerically. The 
full one- loop result, supplemented by the leading two-loop contributions taken over 
from the real MSSM are implemented into the public Fortran code FeynHiggs2.0. 
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1 Introduction 



The search for the lightest Higgs boson is a crucial test of Supersymmetry (SUSY) which 
can be performed with the present and the next generation of accelerators. The prediction 
of a relatively light Higgs boson is common to all supersymmetric models whose couplings 
remain in the perturbative regime up to a very high energy scale A precise prediction 
for the mass of the lightest Higgs boson and its couplings to other particles in terms of the 
relevant SUSY parameters is necessary in order to determine the discovery and exclusion 
potential of LEP2 and the upgraded Tevatron, and for physics at the LHC and future 
linear colliders, where eventually a high-precision measurement of the properties of the 
Higgs boson might be possible @. 

The case of the Higgs sector in the CP-conserving MSSM (rMSSM) has been tack- 
led up to the two-loop level by different methods such as the Effective Potential (EP) 
method HH^H, the renormalization group (RG) improved one-loop EP approach f?j] and 
the Feynman-diagrammatic (FD) method using the on-shell renormalization scheme |||9| . 
This method has provided the only complete calculation at the one-loop level includ- 
ing momentum dependence pll and furthermore the dominant logarithmic and non- 



logarithmic corrections at the two-loop level H,ll|. The application of different methods 



lead to thorough comparisons between the different approaches. 

In the case of the MSSM with complex parameters (cMSSM) the higher order cor- 
rections have been performed, after the first more general investigations [I2|, in the EP 



approach and in the RG improved one- loop EP method [ 13 , |lj , |15 , 16| , 17 1 . In the context of 
the FD approach, so far an exploratory study, including a calculation of the leading terms 
has been performed QlS | . A full one-loop calculation, including momentum dependence, 



as well as an evaluation of the leading two-loop corrections have been missing so far. 

This paper provides the next step into this direction: We present the full one-loop 
calculation for the Higgs boson masses and mixing angles in the cMSSM. Concerning the 
Higgs boson masses, the full momentum dependence has been taken into account. In this 
paper we present briefly the calculation and discuss the effects on the Higgs boson masses 
that originate from the new terms that extend the available calculations obtained using 
the EP and RG approach. More details about the evaluation as well as a more detailed 



numerical analysis can be found in Refs. [jl9|, [20fl . All results are incorporated into a public 
Fortran code, FeynHiggs2.0. 

The rest of the paper is organized as follows. In Section 2 we review the Higgs 
sector and the scalar quark sector of the cMSSM, providing the relations of physical 
and unphysical parameters, the masses and the mixing angles. Section 3 contains the 
numerical analysis. We conclude with Section 4. 

2 Calculational framework 

2.1 The tree-level structure of the cMSSM Higgs sector 

The cMSSM Higgs potential reads [p2T|: 



V = m\H x Hx + m 2 2 n 2 n 2 - (m\ 2 e ah H\H\ + h.c.) 
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where m\, m 2 , m 2 2 are soft SUSY-breaking terms {m\ 2 can be complex), g, g' are the SU (2) 
and U(l) gauge couplings, and e 12 = — 1. The doublet fields Hi and are decomposed 
in the following way: 
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£ is a phase between the two Higgs doublets. In the Higgs potential it appears only in the 
combination £' = £ + arg(m 2 2 ). From the unphysical parameters in Eq. ([!]) the transition 
to the physical parameters (including the tadpoles) is performed by the substitution (see 
Refs. mm for details): 



vi, v 2 , gi, g 2 , mj, m 2 , 



\m 



12 h 



£' -> e, s w , Mz, tan/?, M|±, T h , T H , Ta . (3) 



e is the electric charge, s 



2 

w 



with cw = M w /M z , where M^/ and are the 



masses of the W and the Z boson, respectively. Furthermore, tan (3 is the ratio of the two 
vacuum expectation values, tan/5 = v 2 /v\ (sp = sin/3 , C/3 = cos/5 ), M^± is the mass of 
the charged Higgs boson, and T x ,x = h,H,A denote the tadpoles of the fields h, H and 
A, see below. 

At the tree-level, the tadpoles have to be zero. In the case of Ta this can only be 
fulfilled if £' = 0. Thus CP- violation is absent at the tree-level. The bilinear part of the 
Higgs potential has to be diagonalized to obtain the mass eigenstates. In the CP-even 
sector this is done with the help of the angle a, and results in the two neutral CP-even 
Higgs bosons h and H. In the CP-odd sector the diagonalization can be performed with 
the angle ft, leading to the CP-odd A boson and the neutral Goldstone boson G. The 
charged Higgs sector is also diagonalized with the angle ft, resulting in the charged Higgs 
bosons H ± and the Goldstone bosons G^. One arrives at the following masses at tree- 
level: 
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M H ± and tan ft are chosen as input parameters. The entries for the Goldstone bosons G 
and are to be understood in the Feynman gauge. Since there is no CP-violation in 
the cMSSM Higgs sector at tree-level, there is no mixing between h and H and the fields 
A and G. 
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2.2 The fermion and gaugino sectors of the cMSSM 

Possibly CP-violating parameters occur in all other SUSY sectors of the cMSSM. Most 
important for Higgs boson phenomenology is the scalar quark sector. The mass matrix 
of two squarks of the same flavor, q L and qn, is given by 

Mq ~ I m a X Mn + m 2 ' [b) 



q 1V1 R 



with 



M 



I = M| + M§ cos 2(3 (If — Q q Sw) 



M' 



I = M±, + M 2 z cos2{3Q q s 2 w (6) 
X q = A q — //{cot (3 , tan/3}, 

where {cot/5 , tan/3} applies for {up, down}-type squarks, respectively. A q and /i can be 
complex. In an isodoublet the SU(2) symmetry enforces that Mq has to be chosen equal 
for both squark types. The Mq, on the other hand can be chosen independently for every 
squark type. In the scalar quark sector of the cMSSM N q + 1 phases are present, one for 
each A q and one for \i. The squark mass eigenstates are obtained by diagonalizing the 
mass matrix and are given by 



2 2 , 1 

m <h, 2 = m <i + 2 



Ml + MlT y/(Ml - Ml) 2 + Aml\X q \ 2 } . (7) 



The masses are independent of the phase of X q . However, the phase of A q affects the 
mass eigenvalues since it changes the absolute value of X q . 

The other possibly complex parameters are Mi and M 2 , which are the soft SUSY- 
breaking parameters in the gaugino sector, and rrig, the gluino mass. These parameters 
are less important for the Higgs sector phenomenology; rrig enters the predictions only at 
the two-loop level. 

2.3 The neutral Higgs boson sector at the one-loop level 

The inverse neutral Higgs boson propagator matrix in the cMSSM at the one-loop level 
is given by 



(Amggs) 1 = -i 



±Hh{q 2 ) q 2 -m] I + t HH {q 2 ) ± HA (q 2 ) 

{ t Ah {q 2 ) W<7 2 ) q 2 -m\ + t AA {q 2 ) j 



E denotes the renormalized Higgs boson self-energies (at the one-loop level). CP- violation 
occurs, i.e. mixing between the CP-even Higgs bosons h, H, and the CP-odd Higgs boson 
A occurs if the self-energies T< A h = ^ha and/or T, A h = T,hA are non-zero. This can 
happen if the complex parameters in the cMSSM possess an imaginary part. Also the 
pure CP-even self-energies, E^, T^hh and T>hH, and the pure CP-odd self-energy, T, AA , 
are numerically affected if some cMSSM parameters are complex. Details about the 
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calculation of the renormalized Higgs boson self-energies and the on-shell renormalization 



procedure can be found in Refs. [18. 19 



The poles of the Higgs boson propagator matrix are the squares of the pole masses of 
the three neutral Higgs bosons including higher-order corrections. The masses are denotes 
by m hl , m h2 , m hz with m hl < m h2 < m hz . 

Besides the Higgs boson masses, also the effects of the higher-order corrected self- 
energies on the Higgs boson couplings to SM gauge bosons and fermions can be evaluated. 
In the limit of q 2 = 0, i.e. S(g 2 ) — > S(0), the transition from the tree-level states h, H, A 
to the mass eigenstates at higher orders, hi,h 2 ,h 3 can be described with the rotation 
matrix U, 

' Mn «12 "13 \ / h \ ( h \ 

U2i u 22 u 23 \ \ H \ = U \ H . (9) 
"31 u 32 u 33 ) \ A ) \ A J 

U includes the dominant corrections (coming from Higgs boson propagators) into the 
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effective couplings [22]. The explicit formulas are given in Refs. [16, IS 



The Higgs boson self-energies have been evaluated by taking into account all sectors of 
the cMSSM, including possibly complex parameters and the full momentum dependence. 
The evaluation has been done with the help of the programs FeynArts |23] (using the 



MSSM model file [B^]) and FormCalc p5| . The results, supplemented by the leading two- 
loop contributions taken over from the rMSSM have been transformed into the 
Fortran code FeynHiggs2.0 p6 |. The code evaluates the Higgs boson masses, the mixing 



matrix and the corresponding corrections to the couplings of the Higgs bosons to SM 
gauge bosons and fermions. FeynHiggs2.0 is available at www.feynhiggs.de . 



3 Numerical examples and discussion 

In this section numerical examples are presented and discussed. They are meant to 
illustrate on the one hand the possible effects of complex phases in the MSSM and on the 
other hand the effects of the newly evaluated terms. For a more detailed phenomenological 
analysis constraints on CP- violating parameters from experimental bounds, e.g. on electric 
dipole moments (EDMs), have to be taken into account [p7 |. However, in our analysis 



below we only take non-zero phases for A t = Ab and M 2 , which are not severely restricted 
from EDM bounds. 

The numerical analysis given below has been performed in the "CPX" scenario [28], 
where the parameters are fixed to 

Msusy = 500 GeV, \A t \ = \A b \ = 1000 GeV, /x = 2000 GeV, M 2 = 500 GeV, 
M H ± = 150 GeV, tan/5 = 5 , (10) 

if not indicated differently. The phases are always defined explicitly. However, our analysis 
is confined to M = 0, since this is the most restricted phase, see e.g. Ref. [^{J and 
references therein. The CPX scenario has been defined in order to maximize the effects 
of complex phases. It should be kept in mind that while the relatively high value of is 
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not realized in minimal models like mSUGRA, mGMSB or mAMSB, it can lead to large 
effects especially in the b/b sector. 

In our numerical analysis we concentrate on the Higgs boson masses derived from the 
one-loop corrections only This is sufficient to show the effects of the newly evaluated 
corrections. A more detailed numerical analysis, including also the effects on the rotation 
matrix and correspondingly to the Higgs boson couplings to SM gauge bosons and fermions 



can be found in Refs. 19 



In Fig. |T| we show the two lighter masses, and rrih 2 , for the parameters in Eq. fllOD 
as a function of <pA t - The dotted line shows the result for the (s)fermion sectoiQ in the 
q 2 = approximation, the dashed line shows the full cMSSM with q 2 = 0, and the full 
line also includes the momentum contributions. The effect of the subleading one-loop 
contributions, i.e. the ones beyond the (s)fermion sector, can give rise to changes in 
and m/j 2 of 0(4 GeV). The further inclusion of the momentum dependence can result in 
a shift of of 0(2 GeV) for large values of <pA t - 




Figure 1: and are shown as a function of <pA t f° r the parameters of Eq. (|10|). 

The dotted line shows the result for the (s)fermion sector in the q 2 = approximation, 
the dashed line shows the full cMSSM with q 2 = 0, and the full line also includes the 
effects of the non-vanishing momentum. 



All neutral Higgs boson masses, m^, rrih 2 , rrih 3 , are shown as a function of M#± in 
Fig. |] for the parameters of Eq. fllPl) . It becomes apparent that especially when two 
mass eigenvalues are close to each other, the newly evaluated terms and the momentum 
dependent contribution can change the results by several GeV. Especially the inclusion 

1 With the (s)fermion sector we denote here the t/t/b/b sector. The effects of the other fermions and 
sfermions are numerically small and stay below 1 GeV (l^,|20) . 
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M H ±[GeV] 




Figure 2: m^, m^ 2 and m/, 3 are shown as a function of M#± for the parameters of 
Eq. ([[I]) and 4>A t = n/2. The line styles are as in Fig. [I]. 




Figure 3: A 32 = — rrih 2 is shown as a function of M#± for the parameters of Eq. (|10|). 
The line styles are as in Fig. |l|. The upper and lower curves correspond to <pA t — 7r /2 and 

0A t =O. 
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of the momentum dependence results here in a larger mass gap of the two lightest masses 
compared to the result where the momentum dependence has been neglected. 

Of particular interest is the mass gap between the two heavier neutral Higgs boson 
masses. For values of M H ± > 200 GeV, m^ 2 and rrih 3 can be very close to each other, 
which makes their experimental resolution at a collider experiment difficult. In Fig. we 
show A 32 = m h3 — m h2 as a function of M H ± in the rMSSM, i.e. for 4>A t = (lower set 
of curves) and for (fi^t — 71 / 2 (upper set). The different line styles are as in Fig. |l]. In 
the chosen scenario A32 is always larger in the case of <pA t — 7T /2 as compared to the 
case where <pA t = 0. The induced difference in A32 can be larger than 50 GeV. The 
impact of the non-(s)fermionic terms can be of 0(5 GeV). It should be kept in mind that 
changing the phase of A t also effects the absolute value of X t . In fact, in Refs. flT9| , |2"0|| 
it is demonstrated that every mass gap that can appear in the cMSSM can (for another 
choice of parameters) also be accommodated in the rMSSM. 

Finally, in Fig. |4] we show the effect of 4>m 2 011 the lightest Higgs boson mass in the 
scenario of Eq. flTU|). The dashed (full) line shows the difference of the full MSSM with 
q 2 = (7^ 0) to the (s)fermion approximation with q 2 = 0. The effect of the non- 



<t>M 2 




Figure 4: The effects on of the different sectors normalized to the (s)fermion sector 
evaluation only (i.e. the (s)fermion sector contribution subtracted) is shown as a function 
of 4>m 2 f° r </v = 4>A t — and the other parameters as in Eq. (|T^). The dashed (full) line 
shows the result for the full cMSSM with q 2 = (^ 0). 



(s)fermionic contribution are in this case of 0(3 GeV), and the momentum dependence 
induces a shift of Oil GeV). The effect of 4>m 2 is ^ 1 GeV. Thus the phase of M 2 has 
a much smaller impact than the phases of the trilinear couplings. In a more detailed 



analysis of the cMSSM parameter space (see Refs. [19|,|20[) it has been found that large 
effects of the phases of the gaugino mass parameters on , rrih 2 , vn^ or on the mixing 
of CP-even and CP-odd states can only occur if two mass eigenvalues are very close to 
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each other. Thus, the effect observed in Ref. |15| that the mixing of the two heavier Higgs 
states depends strongly on the phase of or M\ happens only in a very small part of 
the cMSSM parameter space. 



4 Conclusions 

We have presented a complete one-loop calculation of the Higgs boson masses in the 
MSSM with complex parameters. The calculation has been performed in the Feynman- 
diagrammatic approach, using the on-shell renormalization scheme. Besides the full spec- 
trum of cMSSM particles also the momentum dependence has explicitly been included. 

In the numerical analysis we have investigated the effects of the non-(s)fermionic sec- 
tors as well as the effects of the momentum dependence on the three neutral Higgs boson 
mass eigenvalues, m^, and m^. The analysis has been performed for a set of cMSSM 
parameters that maximizes the effects of the CP-violating phases. In this case we find 
that the corrections of the non-(s)fermionic sector can be of (9(5 GeV), while the mo- 
mentum dependence induces a shift of up to 2 GeV. These effects are more pronounced 
if two of the mass eigenvalues are close to each other. The observed effects on m^, rrih 2 
and rrih 3 of the newly included corrections are thus of the order of several GeV and have 
to be taken into account in phenomenological analyses of the cMSSM Higgs boson sector 
in view of the prospective experimental precision at the next generation of colliders. The 
phase of the gaugino mass parameters, 4>m 2 and <pMn are found to have only a relatively 
small effect, except in a very small regions of the parameter space, where two of the mass 
eigenvalues are very close to each other. 

The results of the full one-loop calculation, supplemented with the dominant and 
subdominant corrections taken over from the rMSSM, have been implemented into the 
Fortran program FeynHiggs2.0. The code can be obtained at www.feynhiggs.de . 
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